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We demonstrate that in the atomic-scale limit the thermal conductance /C of the FPU model and 
its variants strongly deviates from the mesoscopic behavior due to the relevance of contact resistance. 
As a result, atomic chains follow log AC = vlogT, where the power law coefficient v is exactly two 
times larger than the mesoscopic value. We smoothly interconnect the atomic and mesoscopic limits, 
and demonstrate that this turnover behavior takes place in other nonlinear FPU-like models. Our 
results are significant for nanoscale applications, manifesting an atomic thermal conductance with 
temperature scaling superior to the mesoscopic limit. 
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The Fermi-Pasta Ulam (FPU) model [l| and its vari- 
ants provide an ideal test-bed for addressing fundamental 
issues in statistical mechanics such as equipartition of en- 
ergy, the onset of chaos in nonlinear dynamical systems, 
and the validity of macroscopic laws in low dimensional 
systems j2j. In particular, the thermal conductivity of the 
FPU model has been extensively investigated, demon- 
strating a breakdown of the normal-diffusional Fourier's 
law dynamics The underlying question addressed in 
these studies has typically been whether a one dimen- 
sional (ID) chain of oscillators with a specific force-field 
can demonstrate equilibrium and dynamical properties 
characterizing macroscopic objects. Thus, simulations 
were mostly carried out using long chains of 10 3 — 10 4 
beads, targeting the mesoscopic- macroscopic limit. 

Nonetheless, in recent years the thermal conduction 
properties of nanoscale junctions has been of fundamen- 
tal and practical interest Here, the physical setup 
includes a nanoscale object, e.g., a nanotube Q or an 
alkane molecule @, coupled to two (or more) thermal 
contacts. As devices size shrinks, an enormously impor- 
tant question is the applicability of the macroscopic dy- 
namical laws at the nanoscale. Even in the framework of 
classical mechanics, contact effects and the constriction 
geometry may imply on the onset of new scaling rules. 
The interest in nanoscale thermal conduction is driven 
by various challenges. One of the major problems in 
molecular electronics is junction heating, limiting the de- 
vice stability. Understating heat transport in molecular 
systems thus attract considerable attention both experi- 
mentally 0,1,0, Hi and theoretically d E3, ISM E3 • 
Managing thermal transport across interfaces [14| is also 
important in microelectronic devices Q, biophysical ap- 
plicatio ns )la| . and in thermoelectric energy conversion 
devices |16j . 

Motivated by these challenges we focus here on the 
atomic-scale steady-nonequilibrium FPU model and its 
variants. As the effect of the contact (interface) and 
internal nonlinear interactions cannot be trivially sep- 
arated, a complex-new dynamics is revealed at different 
temperature domains for both weak and strong system- 



bath interaction strengths. Specifically, using numer- 
ical simulations we demonstrate that when the con- 
tact resistance controls the dynamics, e.g., for extremely 
short chains or at relatively low temperatures (yet above 
the harmonic-anharmonic transition) , the conductance of 
FPU-type systems follows a K, cx T v power law, where 
v = 1/2 for the /3-FPU model. This behavior stands in 
a sharp contrast to the mesoscopic /3-FPU limit where 
the value v = 1/4 is obtained, as expected from the 
phenomenological Debye theory. We justify our results 
within the effective phonon theory [lj, [HI > and further 
interconnect the atomic and mesoscopic limits, demon- 
strating a smooth turnover of the dynamics with increas- 
ing chain size and temperature. Our simulations also 
pinpoint on the critical temperature where contact resis- 
tance is secondary to the bulk thermal resistance. 

We consider a ID lattice of N atoms whose Hamilto- 
nian reads 
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Here Xi is the displacement from the equilibrium position 
of the i — th particle, xq and xn+i are the fixed bound- 
aries, /3 measures the strength of nonlinear interactions, 
and g s (s = 2,4) are Boolean variables taking the values 
of or 1. In what follows we will consider three models: 
(i) A harmonic model (-H2) with g 2 — 1, <?4 = 0, (ii) a 
quartic model (-H4) with g% = 0, 54 = 1, and (hi) the 0- 
FPU model (Hp) with g 2 = 1, 34=1. The quartic model 
is introduced here to facilitate identifying the dynamics 
of the FPU model at high temperatures. 

The direct way to determine the thermal conductance 
JC of a ID chain is to couple the left and right ends of 
the system with two thermal baths at temperatures Tl 
and Tfl, respectively. In our simulations we use Langevin 
thermostats, with the motivation to simulate real exper- 
iments on the nanoscale, where a small molecule con- 
nects to macroscopic solids [1, H|. I n such systems the 
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contact resistance is an unavoidable issue. We note that 
in mesoscopic scale simulations Nose-Hoover thermostats 
are typically used [H, H3, Hj]- However, these ther- 
mostats do not correctly reproduce the weak-coupling 
regime [3(, which is of interest here as well. The chain's 
particles obey the following equations of motion, 



dH 

dxi 



(iR&i - F R ) S i<N - {iL&i ~ F L ) Sis. 
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where Jl,r is the coupling strength between the sys- 
tem and the L,R bath, and the Gaussian thermal white 
noises obey the fluctuation-dissipation relation (Ub = 1), 

(FlA*)) = 0, {F L , R (t)F L ,R{°)) = ^L, R T L , R S(t). The 
long-time heat flux can be calculated between every two 
sites as the force exerted by the i — th particle on the i + 1 
oscillator, Jj = -(xi[g2$l (x i+ i - Xi) + giP(xi+i - Xi) 3 ]} 
where the average reflects the time and ensemble av- 
erage, performed after the system has reached steady- 
state. The thermal conductance is defined as, K, = 
at Etl 1 Ji/(N - 1), where AT = T L - T R . 
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FIG. 1: Thermal conductance of atomic FPU chains (circles 
or dots), a quartic chain (dashed line) and an harmonic chain 
(full line) as a function of the average temperature T — (Tl + 
Tr)/2. y L = jr = 0.6, Q, 2 = O.Of, /3 = 0.01/4. The main 
plot displays data for N=2. The lower (upper) panel displays 
data for TV = 5 (N = 10). 

In Fig. [TJwe present the thermal conductance of rela- 
tively short chains with N = 2, 5, 10, for the harmonic, 
quartic, and the FPU models in the strong friction limit 
using Q = 0.1, 7 = j L m = 0.6, j3 = 0.0025, and 
Tl = O.8T7?. The following observations can be made: (i) 
Disregarding anharmonic interactions, the conductance 
does not depend on temperature, (ii) the quartic model 
shows a power- law dependence, K ~ T 1 / 2 , (iii) the FPU 
model smoothly interconnects the harmonic and quar- 
tic results, with a crossover temperature around 10~ 2 , 
and (iv) for long chains the power law coefficients at 
low and high temperatures are distinct. The upper in- 
set (N = 10) demonstrates this effect, to be discussed in 
more details in Fig. [51 We can explain these results as 
follow. For a ID harmonic chain it can be analytically 



shown [2E 
given by 



that the average stationary thermal flux is 
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where tf>(j) = and e~ a = l + z/2- ^/z + z*/4, 

z = S! 2 /7 2 . In the strong coupling limit, z — > 0, and we 
obtain a ~ y^z, resulting in (f)(1) ~ 1/2 for N = 2 and 
(f)(1) ~ 1 for N ^> 1. Thus, the strong coupling harmonic 
conductance reduces to 
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In the opposite weak-coupling limit (z — > 00) we get a ~ 
lnz, 0(1) ~ 1/z, and K, = 7/2. irrespective of chain size. 
In the present case, Eq. Q predicts the values 0.0041, 
0.0064, 0.0070 for the N = 2,5, 10 chains (respectively). 
Our simulation data perfectly agree with these numbers. 
At low temperatures, since the dominant contribution to 
the transport comes from the quadratic potential, the 
FPU model obeys (TJJ, alike a pure harmonic system. 

We turn next to the low temperature quartic model. 
Here a power- law dynamics is observed with IC ~ T° D 
(N=2). In the high temperature regime (T > 10~ 2 ) 
both FPU and quartic models yield IC ~ T aA . These 
observations stand in a sharp contrast to other results 
observing a T 1 / 4 dependency [13, [2y, [2l|] . Note however 
that these works were concerned with the thermal con- 
ductivity of periodic or long anharmonic chains, where 
contact effects were negligible or non-existing. In con- 
trast, in short atomic systems contact resistance plays a 
crucial role 0, H, We consider its influence next. 

Assuming that the total thermal resistance of an 
atomic-scale chain is given by the sum Rt = Rm + Rc, 
where Rm is the bulk molecular resistance, and Rc is the 
contact resistance, one finds that the total conductance 
satisfies K-t = ICmICc/{IC-m + fcc), where K, n = 1/R n ; 
n = T, M, C . Since the molecular resistance increases 
with size (not necessarily linearly), we can generally as- 
sume that for short chains K~t ~ K-c, while for long sys- 
tems JCt ~ A"Af • Next we estimate separately the contact 
conductance JCc and the molecular contribution Km- 

It is well justified that a sort of a Virial theorem holds 
for nonlinear potentials, resulting in normal mode spectra 
which are pseudo-harmonic for both the quartic potential 
and the FPU model [17j, with the dispersion relation 
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Here the integer k counts the normal modes, ujk 
2 sin(fc7r//V), and the renormalization factor satisfies 



a = 92 + 94 



2/3 (U 4 )g 
« 2 (V 2 )h' 
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The average over the potential energy is taken with re- 
spect to either the H 2 , H 4 or the Hp Hamiltonian. It 
is remarkable that the factor a does not depend on the 
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wavenumber k. The above relation is valid for systems at 
equilibrium. In our simulations only a small temperature 
bias is applied, thus we pose the ansatz that ([5|) holds in 
the nonequilibrium regime as well, i.e., the equipartition 
relation, T = a^i 2 Lof,{q\)T 1 is valid close to equilibrium, 
where averages are taken with respect to the mean tem- 
perature T — (Tjr + Tr)/2. Here qk's are the normal 
modes of a ID harmonic chain. Our numerical results, 
validating flSJ below, justify this conjecture [24[. Calcu- 
lating the thermal averages in ([B]), in the high tempera- 
ture regime the dynamics of the FPU model follows the 
behavior of the quartic model, while at low T it coincides 
with the harmonic dynamics, 
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Since for short chains the boundary resistance dominates 
heat transfer, the conductance is assumed to follow Eq. 
(H|), Kc oc att 2 /j, corrected by the renormalization fac- 
tor ([7]). This leads to the following scaling law 



Kc oc a 



T2 HighT 
T° LowT. 
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This behavior is significant since a a macroscopic- 
phenomenological theory provides a very different scal- 
ing law for the FPU model and its variants: Within the 
Debye formula the heat conductivity can be written as 
k = J2k CkV k Tk > where Cfe, Vk and Tfe are the specific heat, 
phonon velocity, and the phonon relaxation time respec- 
tively. Based on this expression, the following relation 
can be derived [l9l l2fj|. k oc yfa/e, where a is the same 
factor as in Eq. (J7J), and e = (V^n/igiVi + g^V-im- 
In the high temperature limit this results in the follow- 
ing proportion, valid for both the quartic and the FPU 
models, K oc ^/a tx Ti. Assuming that the temperature 
gra dient is linear in our system (besides the contact drop 
|2l|). see [22j], we retrieve the bulk conductance 

JC M ocV^~Ti (9) 

We thus conclude that Kc oc T v whereas Km on T u / 2 . 
While in the low temperature limit and for short chains 
we expect Kc to dominate the overall conductance, in the 
opposite limit the bulk resistance dominates, and even- 
tually the scaling K ~ T 1 / 4 emerges. 

We present next data for a longer chain of iV=20, see 
Fig. Focusing on the quartic model (circles), it is evi- 
dent that the low-T and the high-T power-law coefficients 
deviate. Specifically, at low T the quartic model follows 
K ~ T 0A1 (upper panel) while at high T, K ~ T 026 
(lower panel). The high T values are characteristic for 
the FPU model as well (□). Furthermore, we can system- 
atically extract the power-law coefficients as a function 
of chain size at different temperature ranges, see Fig. [3] 
While at low T the power-law behavior is approximately 
fixed, K~T V with v ~ 0.4 - 0.5, at high T a clear tran- 
sition to the mesoscopic result is obtained with a power 



v ~ 0.25 for N > 20 at T > 1. For convenience the slopes 
were estimated using the quartic-potential data. 

We can further identify the critical temperature T c 
where bulk conductance dominates contact effects. We 
conjecture that R T ~ aT" 1 / 2 + where the first 

(second) term reflects the contact (bulk) resistance, and a 
and b are (iV-dependent) constants. A plot of RtT 1 / 4 vs. 
T -1 / 4 should thus become flat above a critical tempera- 
ture. For N — 50, the inset of Fig. [3] demonstrates this 
turnover at the value 1.1, which translates into T c ~ 0.7. 
We estimate that at this value the phononic mean free 
path is significantly shorter than the molecular length. 
The coefficients a and b are roughly the slope of the linear 
line and the asymptotic (constant) value, respectively. 




FIG. 2: Power law behavior at low and high temperatures. 
1L = jr = 0.6, Q. 2 = 0.01, P = 0.01/4, N = 20, an FPU chain 
(□), a quartic chain (o). The upper and lower panels zoom 
on data in the low and high temperature limit respectively. 
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FIG. 3: Power low behavior of the quartic potential as a func- 
tion of chain size in the high temperature regime, T > 1 (full 
line), low temperature regime T < 0.01 (dashed). Data was 
obtained by studying the slope of log-log plots of conductance 
vs. temperature, with parameters as in Figs. [1] and [2] In- 
set: Thermal resistance times T 0,26 for JV=50, marking the 
dominance of bulk conductance for T > 0.7. 

The divergence of the bulk conductance from the con- 
tact conductance, given by the scaling laws ([S]) and © 
for the /3-FPU models, should take place in other mod- 
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els of confining potentials. In particular, we study the 
thermal conductance of the 2-6 potential, H = Ht + 

g 2 Sl 2 V 2 + g4f3V 6 , with V 6 = gEilo^i+i _ x >-) 6 - Here 
Ht is the kinetic energy and V% is the harmonic term, 
same as in ([2"]), In the high temperature limit we can 

estimate a by [17| a ~ = 3r ^ /6) T 2 / 3 , where T 

is the Gamma function. Short chains, controlled by the 
contact resistance, are expected to follow K. oc a ~ T 2 / 3 , 
beyond the harmonic-anharmonic transition. In contrast, 
long chains should obey K. oc \fo. ~ T 1 / 3 . Our numerical 
simulations have verified this behavior (not shown) . Sim- 
ilarly, studying the a-FPU model, we obtained a power 
law behavior with v = 1/3 for short chains, contrasting 
the mesoscopic v — 1/6 value [201 ] . 
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FIG. 4: Conductance of the harmonic chain (full), quar- 
tic chain (dashed) and the /3-FPU chain (dotted) for vari- 
ous chains in the weak coupling limit 7 < fi, Q, 2 — 0.01, 
Jt = 7fl = 0.05, /3 = 0.01/4. 

As our scaling argument relays on the validity of Eq. 
([3|) for nonlinear systems, renormalizing the phonon spec- 
trum, we further test the applicability of this relation 
in the weak coupling limit, 7 < SI. In this parameter 



range this expression predicts that the conductance of 
both harmonic and FPU- like systems should be constant, 
independent of temperature and size. Our numerical sim- 
ulations, Fig. [31 indeed confirm this expectation: At 
low temperatures both the harmonic potential and the 0- 
FPU model yield K ~ 7/2 = 0.0025. In the high temper- 
ature limit small deviations from this value disclose con- 
tributions to energy transfer beyond the effective-phonon 
theory [2l[. It is also worth noting the intricate dynam- 
ics for 7V=20. As expected, with increasing tempera- 
ture the FPU model nicely interpolates the harmonic to 
the quartic limit. However, quite interestingly, since the 
conductance of the anharmonic chain lies below the har- 
monic value, the FPU model demonstrates a decrease of 
JC with increasing T around the harmonic-anharmonic 
transition, see the rightmost panel in Fig. ((J). 

In summary, we presented here new results for the tem- 
perature scaling of the conductance of anharmonic chains 
with confining potentials. In the atomistic limit we found 
that the contact thermal resistance controls the junc- 
tion conductance, leading to a dynamics significantly dis- 
tinct from the mesoscopic limit where bulk conductance 
dominates. These results are significant fundamentally, 
demonstrating a new atomic limit with enhanced ther- 
mal conductance for the celebrated FPU model, driven 
by interface effects. From the practical point of view our 
results deliver an encouraging message for nanoscale ap- 
plications, manifesting that the combination of contact 
effects and nonlinearity could be beneficial for thermal 
transport, leading to an enhanced conductivity in com- 
parison to the macroscopic limit. 
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